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ABSTRACT 

We have investigated the dipping activity observed in the high-mass X-ray binary 4U 1907+09 and shown that the source continues to 
pulsate in the "off" state, noting that the transition between the "on" and "off" states may be either dip-like or flare-like. This behavior 
may be explained in the framework of the "gated accretion" scenario proposed to explain the flares in supergiant fast X-ray transients 
(SFXTs). We conclude that 4U 1907+09 might prove to be a missing link between the SFXTs and ordinary accreting pulsars. 

Key words, pulsars: individual: - stars: neutron ~ stars: binaries 



1. Introduction 

The persistently active high-mass X-ray binary system 
4U 1907+09, discovered in the third Uhiirii survey ( ^Giacconi| 
|et al. 1971[l, consists of an X-r ay pulsar with a spin period of 
~ 437.5 ( [Makishima et al.|1984) l and a highly reddened compan- 
ion in an eccentric e ~ 0.28 orbit with a period of P ~ 8.3753 d 
(|in 't Zand et al.|1 998 ). Optical and infrared observations ( |Cox| 
|et al.|2005[|Nespoli et al. 2008) suggest an 08-09 la type super- 
giant donor with L ~ 5 x 10^ Lq and a mass-loss rate of about 
7 X lO^^Moyr The supergiant nature of the companion and 
the strong X-ray variability indicate that accretion most likely 
proceeds from a stellar wind. The distance to the source is esti- 
mated to be between 2 and 6kpc ( Cox et al.|2005 Nespoli et al. 
|2008p , therefore, the X-ray luminosity of the source is uncertain. 
Assuming a compromise distan ce of 4 kpc, the observ ed persis 
tent flux of 

luminosity of about 10-"' 



10 '"erg cm ^ s ' 
ergs"'. 



( in 't Zand et al. 1998 1 implies a 



The broadband X-ray spectrum of 4U 1907+09, similar to 
other accreting pulsars, can be modeled with an absorbed cut- 
off power law. Based on Ginga data, Mihara ("1995 ) reported a 
cyclotron resonance scattering feature (CRSF) at 19keV, which 
was later confirmed with BeppoSAX ( Cusumano et al. 1998| l, who 
also reported a harmonic at ~ 39 keV and a narrow iron emission 
line at 6.4 keV. The CRSF harmonic, however, was not detected 



either in RXTE (i n 't Zand et al.|1998] l and Suzaku ( [Rivers et"ar 
|2010) observations. 

One signature feature of 4U 1907+09 is its dipping activity. 
During the dips, the X-ray flux drops abruptly by a decade or more 
for a period of few minutes to several hours ( |in ' t Zand et al. 1 1 998| 
[Roberts et al. [2001 [[Rivers et al.|2010| |. Similar behavior has been 
reported for other HMXBs, notably Vela X- 1 ( Kreykenbohm et al. 
[2^[Doroshenko et al.[20"TT] i and GX 301-2 ( [Gogiig et al.[201 1^ . 
However, 4U 1907+09 spends a substantial fraction of time (up 
to 60%, ,§ahiner et al.|2012 l in the "off" state, whereas dips are 
rare events in other systems. Unlike Vela X-1 and GX 301-2, no 
pulsations were detected during the short dips of 4U 1907+09 
observed with RXTE ( [in 't Zand et al.|1998[ ). |Roberts et al.[ ( |200T) 
were able to detect pulsed emission during an extended low-flux 
episode observed with AS'CA. Two scenarios have been invoked 
to explain the dips: 1) obscuration of the neutron star by a dense 
wind clump, and 2) cessation of the accretion owing a wind den- 



sity drop. Several authors have suggested that in most cases there 
is no evidence of an increase in the absorption column during 
the di ps ( [in 't Zand et al.|1998[ [Roberts et al.|2001[ [Rivers et al. 
20101, so the former scenario is unlikely. 

In this study, we have examined the nature of the dips in the 
light-curve of 4U 1907+09 in great detail, using recent Suzaku and 
RXTE observations. First, we compared the flaring behavior of the 
source with the one of other HMXBs, and then we attempted to 
link the observed phenomenology with the flaring activity typical 
of the SFXTs. We affirm that similarities of 4U 1907+09 between 
both normal wind accreting pulsars and SFXTs suggest the source 
might constitute a missing link between the two classes. 



2. Observations and facts 

Pulsations during tlie "off-states. For a deeper understand- 
ing of the nature of the source, it is essential to clarify whether 
the residual X-ray emission of the "ofF'-states is pulsed or not. 
This is indeed the case for Vela X- 1 ( [Doroshenko et al.|2011j l 
and GX 301-2 G6gii§ et al. (2011 1. As already mentioned, 
|in 't Zand et al.[ ( [l998[ ) did not detect pulsations during the dips 
of 4U 1907+09 using RXTE observations, however, the source 
lies close to the galactic ridge, and therefore systematic effects 
related to a stronger background might affect observations with 
collimating instruments. On the other hand, ASCA observations 
hint at the pulsations during the two lOks-long extended low-flux 
episodes ( Roberts et al.[[200T i when the source flux was com- 
parable to the "ofF'-states. It is unclear, however, whether the 
pulsations were significant and whether such a low-flux episode 
could be caused by the same mechanism as the short dips. To clar- 
ify this, we analyzed a 123 ks Suzaku observation (ID 401057010 
from May 2006) of the source. We repeated the analysis carried 
out by Rivers et al. (2010") who used this observation to study the 
behavior of the source in the normal flux state, focusing, however, 
exclusively on the source properties during the "off-states". 

Suzaku is equipped with two instruments: the X-ray spectrom- 
eter (XIS), consisting of focusing four units and a collimating 
hard X-ray detector (HXD). To increase the signal-to-noise ra- 
tio, we only used XIS data, combining all four units for timing 
analysis. For data reduction, we used the HEASOFT 6.1 1 anal- 
ysis package, and the set of calibration files v.20111109. The 
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Fig. 1. The observation-long XIS lightcurve (background-subtracted with 4 XIS units combined) and close-up of some of the "off" 
states, all of which are marked with red points in the lightcurve. 
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Fig. 2. Count-rate distribution in the background- subtracted XIS Fig. 3. Background subtracted "on" (top) and "off" (bottom) XIS 
lightcurve (44 s bin size). Two roughly log-normally distributed pulse profiles in 1-10 keV energy range folded with the same 
peaks corresponding to "off" and "on" states can be identified. parameters. 



XIS lightcurve of the entire observation is presented in Fig.[T] 
Several dip episodes are observed throughout the entire observa- 
tion. Similar to Vela X-1 ( [Doroshenko et al.|2011[ ), the observed 
count rate roughly follows a log-normal distribution for both the 
"on" and "off" states, with different mean values as is also evident 
from the lightcurve (see Fig.|2]). Here we define "off"-states as 
the abrupt flux drops to a level of ~ 1 cts/s per XIS unit for at 
least one spin cycle of the pulsar. Note, that it is expected to find 
that the "off"-states contribute to the tail of the "on" flux distribu- 
tion rather than forming a separate peak, should they come from 
random fluctuations in wind density or velocity. 

For the timing analysis, the lightcurve was corrected for the 
orbital motion in the solar and in the binary system using the 
orbital parameters by in 't Zand et al.| ( 1998 1. Phase-coherent 
timing analysis, using the XIS and HXD Ught-curves, yielded a 
constant pulse period of P=441.10(5) s, consistent with the value 
reported by Rivers et al. ( 20 1 1. 

A blind search for significant pulsations in the "off" lightcurve 
is challenging because it only contains a few pulse cycles sep- 
arated by many gaps of various lengths, which alias with the 
intrinsic variability and result in very noisy periodograms if con- 
ventional algorithms are used. On the other hand, folding the "off" 
lightcurve with the period found from "on" lightcurve produces 
a pulse profile with a similar shape (see Fig. |3]l, which is also 
similar to the ASCA folded lightcurve for the extended low-flux 
episode (see Fig. 5 of Roberts et al.|2001 1. To justify the signifi- 
cance of the pulsations one can estimate the probability of getting 
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Fig. 4. The distribution of squared deviations in the synthetic 
pulse profiles from the constant (black) and the best-fit probability 
density function to the distribution (red). Dashed line indicates 
the observed squared deviation for the "off" pulse profile. 



the observed deviation from a constant over the pulse profile 
by chance (for known pulse-period and pulse shape). To do this 
we simulated lO'* modified "off' lightcurves with flux randomly 
distributed with the same average value and dispersion as the 
observed "off" flux. We then folded the simulated lightcurves 
with observed period and calculated the deviation of the resulting 
pulse profile from a constant. These turn out to be distributed 
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log-normally, so one can estimate the probability of getting the 
deviation equal to or larger than observed. We find that it is fairly 
low at Po ~ 6 X 10"^ as summarized in Fig.|4j i.e. the observed 
pulsation is significant. 

"Interdips" or flares It is interesting to observe that the frac- 
tion of time the source spends in the "off" state is so long that 
occasionally two consecutive dips are separated by just a short 
"on" interval of a comparable length to that of the dips. This in- 
terval would be considered a transient flaring episode, should 
only the part of the lightcurve shown in the third upper panel of 
Fig.[T]be observed. We searched for similar events in the source 
light-curves in RXTE archive, that had extensively monitored the 
source. Searching for uninterrupted pointings, which fully cov- 
ered intervals between consecutive dips, we found the two cases 
presented in Fig. [5] We would like to emphasize that the sharp 
transition between the "off" and "on" states produces distinct 
flare-like events in the lightcurve. 

Out of the RXTE archival data (§ahiner et al. 2012) another 
relevant observational fact emerges: the probability of observing 
the source in the "off"-state depends on the orbital phase and is 
significantly lower at periastron. 



3. Discussion 

The discussion section begins with a short summary of the ob- 
served phenomenology: 

- 4U 1907-1-09 exhibits a dipping activity associated with sud- 
den drops in flux by an order of magnitude or more on a 
timescale that is comparable to the spin period. Dips last sev- 
eral minutes to several hours before the flux recovers to the 
"on" value on a similar timescale. 

- 4U 1907+09 is not unique and a similar behavior has been 
observed in other HMXBs, notably Vela X-1 and GX 301-2. 

- In the "ofi^' state, accretion does not cease completely, form- 
ing a distinct peak in the flux distribution while the source 
continues to pulsate. 

- The pulse profile changes significantly during the "ofF'-state. 

- The probability of observing the source in the "off" state 
depends on the orbital phase and is lower at periastron. 

- The fraction of time 4U 1907+09 spends in the "ofl'"-state is 
long enough to make inter-dip intervals occasionally appear 
as flares. 

First, emphasis should be placed on the phenomenological sim- 
ilarity between 4U 1907+09 and the SFXTs. The inter-dips or 
flares shown in Fig|5]exhibit timescales and amplitudes similar 
to the smaller flares observed from the SFXTs. Indeed, we have 
measured an average unabsorbed flux of ~ 4.5 x 10"" erg cm^^s"' 
and ~ 4.5 x 10 '"ergcm^^s"' for the "off" and "on" fluxes re- 
spectively, implying ^ 8.6 x 10^"^"^^ erg s ' for an assumed 
distance of 4 kpc, so is comparable to values observed during the 
flaring periods in SFXTs ( [Sguera et al.|2008[|Bozzo et al.|2008| l. 

On the other hand, 4U 1907+09 in "off" state is brighter than 
the SFXTs in deep quiescence, and it has shorter "off" duty cycle. 
The latter depends, however, on the orbital phase, hence on the 
distance from the companion. One could anticipate, therefore, 
that the source would remain "off" longer should the system be 
wider On the other hand, the duty cycle of SFXTs is also known 
to depend on the orbital phase, and some of them leave the deep 
quiescence and only start to flare in the vicinity of the perias- 
tron ( Sidoli|2011 1. Making such a system more compact would 
extend the flaring period and prevent if from going to deep quies- 
cence, thus making it similar to the 4U 1907+09 with its relatively 



bright "off" states and lower dynamical range than the prototyp- 
ical SFXTs. Luminosity-related changes of the pulse profiles 
similar to ones in 4U 1907+09 have also been observed in SFXT, 
at least in one case (Romano et al.|2009| l. 4U 1907+09 seems to 
exhibit, therefore, properties that are intermediate between the 
SFXTs and the persistent systems. 

The dipping activity in 4U 1907+09 is not unique among 
wind-accreting pulsars. As discussed by [Kreykenbolun et"ar 
( 2008 1 for Vela X- 1 , the drop in X-ray flux must be intrinsic 
rather than due to eclipses by a wind clump obscuring the neutron 
star. Since no "off" flux was known to emerge from Vela X-1 



and 4U 1907+09, Kreykenbohm et al. (2008 ) and gahiner et al. 
(2012 ) have suggested that the centrifugal inhibition might be re- 
sponsible for the cessation of emission. However, residual pulsed 
flux has been detected in Vela X-1 ( jPoroshenko et aL]|2011| ), 
suggesting that a change in the accretion regime might be re- 
sponsible for the flux drop. Following Burnard et al. ( 1983| l 
and |Bozzo et al. ( 2008| l, Doroshenko et al. ( 2011 1 proposed 
for Vela X-1 that, while normally the accretion proceeds via 
Rayleigh-Taylor instability (RTI), this switches off during the 
dips, and Kelvin-Helmholtz instabilities (KHI) become relevant. 
Note that the switch between two accretion regimes naturally ex- 
plains the observed bimodal flux distribution. The same scenario 
may also apply to 4U 1907+09. 

It is interesting to discuss why just a few X-ray pulsars do 
switch "off". For SFXTs Bozzo et al. (2008 ) argue that a strong 
magnetic field is required to explain the observed luminosity 
swings and the characteristic accretion rates in different accretion 
regimes. They estimate that RTIs become inefficient for mass-loss 
rates of the optical companion below: 



M_6 < 280P:?a^ 



6 i ZOUZ-^jC/iorfJ^S^MlO 



[1 + 16/?„io/(5/?Mio)] 



-3/2 



(1) 



Here, Rcw and Rmw are the capture and magnetosphere radius, 
respectively, in units of 10"^ cm; any - P^i^^M^j^ is function of 
the orbital period and of the total mass of the system; is the 
relative velocity between the wind and the neutron star in units of 
10^ cm s"' ; and P53 is the spin period in units of 1000 s. Assuming 
for 4U 1907+09 the parameters reported by |Cox et ar] ( |2005^ , we 
obtain an upper limit for the magnetosphere size of Rm\o < 0.6. 
Routinely observed regime transitions imply, however, that the 
^Mio is close to this value, i.e. large as well. 



Bozzo et al. ( 2008| also provide an estimate for the leak 



rate and the corresponding X-ray luminosity if plasma enters the 
magnetosphere via KHI. For 4U 1907+09 the estimate is 

Lkh - GMns^kh/^ns - 



1 +Pi/Pi 



ergs 



Where pi e are the densities within and outside of the magneto- 
sphere. According to [Bozzo et al.| ( [2008j ), t/kh ~ 0.1, and the 
density ratio is estimated to be between 



^^pJp'c _ (03i]KHh-^RlilP:. 



1 +Pi/Pe 



0.1?7kh/!"'^mio^8 



where h is the fractional height of the area where the plasma and 
the magnetic field coexist, in units of the total thickness of the KHI 
unstable layer ( [Burnard et al.|1983| l. In the case of 4U 1907+09, 
for the observed "off-state" luminosity of ~ 8 x lO^'^ergs \ 
Rm\o ^ 0.25 can be estimated if the KHI unstable layer is rela- 
tively thin {h ~ 0.05), or Rmio > 0.55, if h ~ I, as suggested by 
[Burnard et al.[ ( [T983] l. The latter is in good agreement with the 
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value deduced from stability criteria. For the observed X-ray lu- 
minosity, this implies a magnetized neutron star with B ~ 10^^ G. 

[Kreykenbohm et al. (2008 ) proposed that the magnetosphere 
might grow enough to inhibit the accretion even for moderately 
magnetized neutron star if the wind density in the vicinity of 
the neutron star drops by factor of ~ 1000. This would allow 
reconciliation of the "off"-states with a magnetic field intensity 
of B ^ 2 X 10'^ G like the one deduced from the observed CRSF 
energy. The winds of young supergiants are known to be struc- 
tured, and hydrodynamical simulati ons have predicted a wind 
density contrast as high as 10^ (Os kinova et al .112007 1. These 
predictions are not confirmed by observations, however, because 
the observed X-ray flux typically varies at most by a factor of 
10. Moreover, the "o ff-state" duty cycle o f 4U 1907+09 depends 
on the orbital phase ( [gahiner et al.|201 2^; i.e., the fluctuations in 
wind density which trigger the "off"-states must be comparable 
to the ones associated with the orbital motion. The average wind 
density in 4 U 1907+09 is expected in the CAK model ( [Castor| 
et al.|1975 1 to vary just by a factor of ~ 5 along the orbit, which 
is consistent with the observed orbital lightcurve. The "off"-states 
in 4U 1907+09 are therefore most likely triggered by minor fluc- 
tuations in wind density, so a large magnetosphere would imply a 
magnetic field for the neutron star of B ~ 10'^ G. Such a strong 
magnetic field might help also in explaining the long spin-period 
of 4U 1907+09, although its comphcated evolution definitively 
deserves a dedicated study. 

The observed CRSF energy in 4U 1907+09 implies, how- 
ever, an order of magnitude weaker field B ~ 1.5 x 10'^ G 
in the line-forming region. Considering that the details of the 
magnetosphere-plasma interaction are far from being fully under- 
stood, and that there is significant uncertainty in the distance to 
the source (which reflects in the estimate of the accretion rate and 
consequently of the magnetosphere radius), it may well be that the 
scenario outlined above will also hold for a weaker field. However, 
4U 1907+09 is a very unusual HMXB, and we tend to believe that 
there is something distinctive behind its dipping behavior, and 
the strong magnetic field of the neutron star is a natural candidate. 
As previously discussed for the case of Vela X-1 and GX 301-2, 
the CRSF energy provides an estimate for the magnetic field in a 
line-forming region, and if the line forms far above surface of the 
neutron star, it might underestimate the surface field by orders 
of magnitude. The dipole component of the field changes with 
~ R^^, so scattering in the upper parts of the accretion column, 
or in the accretion stream at heights comparable to the neutron 
star radius atop the polar caps may help explain the observed 
CRSF energies even for magnetar-like surface fields. Note, that 
the complicated observed shapes of the pulse profiles in bright 
accreting pulsars cannot be explained with the direct emission 
from compact polar caps (Kraus et al._2003) , and this possibility 
should be considered irrespectively of the surface field. 

In 4U 1907+09, the luminosity is probably too low to form 
an extended accretion column ILyubarsky & Sunyaev| [T988] l, 
however, scattering of X-ray photons from the polar caps in the 
upper accretion stream may be important. This was discussed by 
|Kraus et ar] ( |2003| l, who find that scattering in the upper accretion 
flow might afi'ect the formation of the pulse-profile and spectra 
in a major way and, in fact, be responsible for the majority of 
the hard photons coming from the source. The magnetic field 
in scattering region might already be compatible with the one 
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Fig. 5. Two SFXT-like flares observed by RXTE and centered 
respectively at ~MJD 55 169.5 and ~MJD 51074.49. In the latter 
case the source remains "ofi"' after a gap in the data for ~ 3 h. 

with respect to the observer will also affect line formation because 
only part of the emission or scattering region will be visible. 

Another scenario that could explain the difference between 
SFXTs and normal HMXBs, invokes properties of the stellar 
wind. Recently, 'Ikhsanov & Finger ( 2012| l have suggested that 
the picture of interaction of the accretion flow and the magneto- 
sphere outlined by |Bozzo et al.| p008| l should be altered if the 
wind plasma itself is magnetized, and this might help explain 
the long pulse periods in some pulsars. We note here that the 
wind magnetic field might also affect the stability of the accretion, 
although how still needs to be clarified. In this scenario the strong 
magnetic field of the optical companion could, in principle, be 
responsible for the difference between the ordinary HMXBs, the 
intermediate systems like 4U 1907+09, and the SFXTs. Such a 
scenario requires further theoretical investigation, while system- 
atic survey of the magnetic fields of the primaries of different 
source classes may be carried out already now to verify whether 
they differ significantly. 



estimated from the observed CRSF energy. Nishimura ( 2008 1 
consider the CRSF formation in an accretion column with a 
height of several kilometers and suggest that line-like features 
are still present in the spectrum. As we discussed for GX 301-2 
( [Doroshenko et al.|2010] l, the orientation of the emission region 
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